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ABSTRACT 

Wc have obtained the most extensive and most accurate photometric data of a Blazhko 
variable MW Lyr during the 2006-2007 observing seasons. The data within each 0.05 
phase bin of the modulation period {Pm = fm^) cover the entire light cycle of the 
primary pulsation period (Pq = fo )■, making possible a very rigorous and complete 
analysis. The modulation period is found to be 16.5462 d, which is about half of that 
was reported earlier from visual observations. Previously unknown features of the mod- 
ulation have been detected. Besides the main modulation frequency /m, sidelobe mod- 
ulation frequencies around the pulsation frequency and its harmonics appear at ±2/ni, 
±4/i„, and ±12.5/ni separations as well. Residual signals in the prewhitened light curve 
larger than the observational noise appear at the minimum-rising branch-maximum 
phase of the pulsation, which most probably arise from some stochastic/chaotic be- 
haviour of the pulsation/modulation. The Fourier parameters of the mean light curve 
differ significantly from the averages of the Fourier parameters of the observed light 
curves in the different phases of the Blazhko cycle. Consequently, the mean light curve 
of MW Lyrae never matches its actual light variation. The $21, ^31 phase differences 
in different phases of the modulation show unexpected stability during the Blazhko 
cycle. A new phenomenological description of the light curve variation is defined that 
separates the amplitude and phase (period) modulations utilising the phase coherency 
of the lower order Fourier phases. 

Key words: stars: horizontal branch - stars: variables: other - stars: individual: 
MW Lyr - stars: oscillations (including pulsations) - methods: data analysis - tech- 
niques: photometric 



1 INTRODUCTION 

The Blazhko modulation of RR Lyrae stars, a phenomenon 
known for about a century, is still one of the open questions 
in astrophysics. The problems and the inconsistencies of the 
existing models wi th observationa l facts were recently dis- 
cussed in details bv lStothersI (|2006r ). therefore we only briefly 
mention here that none of the suggested models (magnetic 
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oblique rotator, IShibahashil (l200d); resonant exc itation of 
nonradial modes, iDziembowski fc Mizerskil (I2OOJ)) can ex- 
plain the complexity of all the observed properties of the 
modulation. 

Alternatively, IStothersI (|2006l ) suggests that the modu- 
lation may be explained by continuous amplitude and pe- 
riod changes of the pulsation due to the action of a turbu- 
lent convective dynamo in the lower envelope of the star, 
with dynamo cycle identical with the Blazhko period. This 
explanation does not involve any nonradial mode compo- 
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nent, the phenomenon is interpreted in the framework that 
Blazhko RRab sta rs are pu r ely fu ndamental mode radial 
pulsators. However. IStotherd l)2006l ) gives only a qualitative, 
rough picture of his model, that have to be checked both 
observationally and theoretically, thus the investigation of 
the phenomenon still remains an important and valid task. 

A lot of efforts have been already made in studying the 
light curve changes of RR Lyrae stars but most of the avail- 
able photometries of Blazhko variables have some defects: in- 
accuracy (visual, photographic data), biased data sampling 
(the photoelectric observations focused mostly on the rising 
branch, maximum phase of the light curve), data inhomo- 
geneity due to rare data sampling on a too long timebase 
(the time scale of period changes and changes in the modu- 
lation properties can be as short as a few years), etc. Most 
of the recent CCD and/or photoelectric observations of indi- 
vidual Blazhko stars are not extended enough to study the 
modulation properties in full detail. The first multicolour 
photometric data that covered each phase of both the mod- 
ulation and the pulsation within a short time-base (o ne sea- 
son) were published for RR Gem (| Jurcsik et al.ll2005l ). Dur- 
ing the period of the CCD observations RR Gem showed, 
however, small amplitude modulation, that limits the avail- 
able information due to the low S/N ratios of the modulation 
signals. 

We started a systematic search for previously unknown 
Blazhko variables at Konkoly Observatory in 2004 using the 
advantage o f full access to an automated 60 cm telescope 
llS6doill2007). The aim has also b een to clear up questionable 
cases. In lSodor fc JurcsikI (l2005f) we ha ve revised the list of 



known Blazhko variables (|Smithlll995f ). The modulation of 
MW Lyrae w as found to be ambiguous as the photographic 
observations (|Gessneilll966 ) did no t confi rm the modulation 
detected in visual data by iMandell (|l970l ). Our CCD obser- 
vations of MW Lyrae affirm the light curve modulation of 
the star , but with a modulation period about half of that 
iMandell (|l970l ') announced. 

Spectroscopic observation of MW Lyr has never been 
obtained, however, based on its short period {Ppuie = 
0.397 days) it should be a relatively hot and metal rich RR 
Lyrae star. 

In the 2006 and 2007 seasons we obtained 
CCD photometric observations of MW Lyrae 
(a = 18'n9'"53'!8,5 = -h31°58'54", J2000). This is the 
first multicolour photometric data set of a large modulation 
amplitude Blazhko variable that is condensed, extended, 
and accurate enough to detect previously unknown 
properties of the modulation. 

In this paper the photometric data of MW Lyr are 
published and analysed. The light curve solution (frequency 
analysis) is discussed using mostly the V band data. A sec- 
ond paper is going to be devoted to the study of the colour 
behaviour during the course of the modulation cycle. 



2 DATA 

The major part of the observations were obtained with the 
automated 60 cm telescope of the Konkoly Observatory, 
Svabhegy, Budapest equipped with a Wright Instruments 
750 X 1100 CCD camera (FoV 17' x 24'). Measurements 
were taken on 177 nights between May 2006 and Jul 2007. 
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Figure 1. Map of MW Lyr and the comparison (CO) and check 
stars (C1,C2, and C3). 



CCD observations with the 60 cm telescope of the Michi- 
gan State University equipped with an Apogee Ap47p CCD 
camera (FoV 10' x 10') and with the Im RCC telescope of 
the Konkoly Observatory equipped with a Princeton Instru- 
ments VersArray 1300B CCD camera (FoV 6' x 6') were ob- 
tained on 6 and 3 additional nights in July 2006 and in Aug 
2007, respectively. Johnson-Cousins BVIc filters were used 
in all the observations. Altogether 5700-5800 data points 
in the BVIc passbands were gathered. Observations in Re 
band were also obtained on 15/3 nights with the Konkoly 
60/100 cm telescopes. 

Data reduction was performed using standard IRAF[j 
packages. Aperture photometry of MW Lyr and several 
neighbouring stars were carried out. In the analysis mag- 
nitudes of MW Lyr relative to GSC2.2 N0223233663 (CO) 
are used. Table [l] lists the basic data and Fig. [l] shows the 
positions of our comparison and check stars. The magni- 
tude differences of the comparison and check stars remained 
constant within 0.010-0.015, 0.009-0.013, 0.008-0.012, 0.009- 
0.013 mag in the B,V,Rc, and Ic bands, respectively (see 
also Fig. [7]). 

Transformation to the standard system was done using 
the B, V, Re, and Ic magnitudes of the stars in the field of 
MW Lyr observed by A. Henden with the USNO Flagstaflt 
Station 1.0 m telescope equipped with a SITe/Tektronix 
1024 x 1024 CCD. A complete list of the positions and stan- 
dard BV{RI)c magnitudes of these stars are available on- 
line as Supplementary Material (Table la). To calculate the 
colour terms of the transformations the nightly instrumental 
V and B light curves were fitted by different order Fourier 
sums to determine the V magnitudes at the moments of the 
B,Rc,Ic observations and the B magnitudes at the mo- 
ments of the V observations. No transformation was applied 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Data of comparison and check stars 



Name 


GSC2.2 ID 


USNO-Bl.O ID 


coord 


y[mag]* 


B 


- y[magr 


V 


-iJcImag]* 


V 


-/c[mag]* 








RA 


Dec 
















CO 


N0223233663 


1219-0328347 


18 19 56.80 


-1-31 58 13.3 


11^841 




0^582 




0^315 




a" 666 


CI 


N0223233659 


1219-0328429 


18 20 06.56 


-1-31 58 56.4 


12^986 




0™682 




0^383 




a" 788 


C2 


N0223233671 


1219-0328258 


18 19 48.61 


-1-31 57 14.6 


13™921 




0^858 




0™480 




ffi"950 


C3 


N02232332613 


1219-0328267 


18 19 48.91 


-1-31 56 54.4 


14™ 363 




0^603 




0^338 




ffi"671 



Standard magnitudes measured by A. Hendcn 



Table 2. Time series of the V magnitude and colour differences 
of MW Lyrae relative to the comparison star CO 



HJD-2 400 000 


AV 


A(B-y) 


A{y-iJc) 


A{y-/c) 


obs* 


53887.34009 


2.218 


-0.075 


-0.064 


-0.041 


1 


53887.35249 


2.227 


-0.073 


-0.035 


-0.035 


1 


53887.35867 


2.250 


-0.087 


-0.020 


-0.020 


1 


53887.36489 


2.257 


-0.085 


-0.020 


-0.023 


1 



* 1) Konkoly 60cm; 2) MSU 60cm; 3) Konkoly Im 

Table 3. Maximum timings and brightness values derived from 
the V light curve. 

maximum time V maximum brightness 
HJD-2 400 000 relative magnitude to CO 



53887.520 
53901.447 
53903.428 



1.501 
1.285 
1.480 



in the /c band of the Konkoly 60 cm data as no colour de- 
pendency of the differences between the instrumental and 
standard magnitudes was found in this band. 

Taking into account the proximity of the comparison 
star to the variable only second order extinction correction 
were applied in the B band. 

B — V, V — Re, and V — Ic colours were derived utilizing 
the V observations and fitted values of the B, Re, Ic curves 
according to nightly Fourier fits for the moments of the V 
measurements. Differential V magnitudes and B — V , V — 
Re, and V — Ic colours of MW Lyr with respect to CO are 
given in TableU] The HeUocentric Julian Date, AV, A{B - 
V), A{V — Re), and A{V — Ic) relative magnitudes, and 
the observatory ID are given in Column 1-6, respectively. 
The entire list of photometric data is available online as 
Supplementary Material. In Table 2a, 2b, and 2c (electronic 
only) the AS, A_Rc and A7c time series are given. 

Maximum timings and maximum brightness values of 
the V data set were determined for 88 epochs (Table [S}. 
The complete list of these data is available online. 



3 RESULTS 

The V light curve of MW Lyrae and the data folded with 
the pulsation and modulation periods are shown in Fig. [S] 
The light curve is strongly modulated, the full amplitudes of 
the amplitude and phase modulations are larger than 0.45 
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Figure 2. Delta V magnitudes versus Julian Date (panel A), 
and data phased with the 0.397674 d pulsation (panel B), and 
the 16.546 d modulation periods (panel C) are shown. 




-0.25 0.25 

pulsation phase 

Figure 3. Residual V light curve of MW Lyr after removing the 
pulsation components of the light curve solution from the data. 
The residuals have larger amplitude at around the minimum, ris- 
ing branch phases of the pulsation than around maximum bright- 
ness due to the large amplitude of the phase modulation. 
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Table 4. Light curve solutions 



Solution 



No. of frequencies Frequency identification 



A (triplet) 

B (quintuplet) 

C (septuplet) 

D 

E 



36 
53 
58 
66 
96 



fc/o,fc : 



12;fc/o-/^,fc = l, 



0; fc/o + /m 


fc=l,...13;/ir 










2/m,fc = l, 


...7;fc/o + 2/m 


,fc = 


1, 


...9 


2/m 


4/m,fc = l, 


2,3;fc/o + 4/n, 


,fc = 


= 1, 


2 




12.5/m,fc = 


2,3;fc/o + 12. 


5/m 


k 


= 1, 


2;fc/o 



frequencies in solution A and fc/o 
frequencies in solution B and fc/o 
frequencies in solution C and fc/o 
frequencies in solution D and 30 additional frequencies (see details in the text) 



/ml '^ — l; 2, 3; /^ 



/o = 2.514621 cd-^ /m = 0.060437 cd-^ //„ = 0.00197 cd" 



mag and 0.07 phase of the pulsation (~ 40 min), respectively 
(see also Fig. [9]). The modulation seems to have larger am- 
plitude at around minimum and rising branch phases than 
at maximum light as Fig.|3]shows. In this figure the residual 
light curve is phased with the pulsation period after the pul- 
sation components are removed. The large amplitude of the 
residuals around minimum phase is the consequence of the 
large amplitude of the phase modulation. In reality, the am- 
plitude of the maximum brightness variation is larger than 
the amplitude of the minimum brightness variation as the 
top and bottom envelope curves of panel C in Fig. [2] show. 

The elements of the pulsation and modulation are: 
TmaxpuiB = 2 453 963.4950 [HJD] + 0.3976742 ■ £puie, 

and 

TmaxBi = 2 453 963.4950 [HJD] + 16.5462 ■ Em- 

The pulsation and modulation periods are those that 
yield the best fit to the V light curve with the pulsation 
(fc/o) and modulation (fc/o ± /m, fc/o ± 2/m, /m, and 2/m) 
frequency components using 'locked' frequency solution (i.e., 
the modulation components are at the positions of the linear 
combination frequencies). Details of the frequency compo- 
nent determinations are given in Sect 3.1. 

Data analysis was performed using the different appli- 
cations of the MUFRAN /TIFRAN packages (JKollathlll990l : 
iKoUath fc Csubrvl 120061 ) . a linear combination fitting pro- 
gram developed by A. Sodor, and the linear and nonlinear 
curve fitting abilities of gnuplolo. 



3.1 The light curve solution 

The Fourier spectrum of Blazhko RRab stars is charac- 
terized by equidistant frequency triplets with frequency 
sep aration identical with the modulation frequency (see 
e.g. ISmith et al.lll999l : I Jurcsik et al]|2005l : iKolenberg et all 
120061 ). In the residual spectr um of RV UMa, e quidistant 
quintuplets were identified bv iHurta et al.l (|2008h . The ap- 
pearance of the modulation frequency (/m) itself was a 
matter of debate for long, but in the spectra of extended 
and accurate datasets fm also shows up unquestionably 
iJurcsik et al.ll2005l . [20061 : iHurta et al.ll2008l) . 

The Fourier spectrum of MW Lyrae is also dominated 
by these frequency components. Pulsation frequencies ap- 
pear up to the 12th order, while the fc/o -f /m and fc/o — /m 
components are present with fc Sj 13 and fc ^ 10, respec- 
tively. The modulation frequency (/m) has an amplitude of 
0.014 mag in the V data, which is as high as the amplitude 
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of the 6th harmonic component of the pulsation. Prewhiten- 
ing the data with the frequencies of the triplet solution sig- 
nificant peaks in the residual spectrum appear, indicating 
that the light curve cannot be accurately fitted simply with 
equidistant frequency triplets. 

Fig. U shows the residual spectra of the V data in the 
— 20 cd~^ frequency range and in the vicinity of the fc/o 
(fc = 1, ...7) pulsation frequencies after prewhitening with 
different frequency solutions. Table 3] summarizes the fre- 
quency solutions applied as more and more frequency com- 
ponents are identified. All the identified frequency compo- 
nents (pulsation plus modulation frequencies) are simultane- 
ously fitted to the original data and this light curve solution 
is used in the next step of the analysis. In Sect 3.3 it is 
documented that the simultaneous fit of all the frequency 
components gives better light curve solution than if the fre- 
quency components are successively fitted and are removed 
from the data in consecutive steps. 

The left panels of Fig. |4]show the residual spectrum af- 
ter removing the triplet frequency solution (solution A). The 
highest peaks appear at fc/o-|-2/m frequencies. The fc/o — 2/m 
and the 2/m frequencies are also detected but with relatively 
small amplitudes. There are altogether 16 frequencies identi- 
fied to belong to the fc/o ± 2/m modulation series. Fitting the 
data with the quintuplet frequencies (solution B) the resid- 
ual spectrum (middle panels in Fig. |3| is still not flat, be- 
sides other peaks two sets of frequencies are evident, one at 
4/m and the other at 12. 5/m frequency separations. Unfortu- 
nately, the frequencies of these series are very close to the ±1 
cd~^ alias components of each other (see their frequency val- 
ues in Table [5]) , that makes the determination of the ampli- 
tudes of these components ambiguous. However, either the 
fc/o±4/m or the fc/o±12.5/m components are removed, mem- 
bers of the other frequency series remain, consequently the 
components of the fc/o ± 4/m and the fc/o ± 12. 5/m frequency 
series are independent signals. In the right panels in Fig. |4] 
it is shown that even if the septuplet solution is removed 
with 5 components of the fc/o ± 4/m series (solution C), sig- 
nals at 12. 5/m separations are still present in the residual. 
These are the highest peaks detected in this spectrum. It is 
not at all clear whether these modulation components are 
indeed connected to the main modulation frequency or it 
happens just by chance that this secondary modulation has 
a modulation frequency very close to 12. 5/m. We can only 
say that the frequencies of this modulation series are within 
the uncertainties at the positions of fc/o ± 12. 5/m. If the 
frequencies of these modulation components are not locked 
to the fc/o ± 12. 5/m positions in the fitting process their 
displacements do not exceed significantly and/or systemat- 
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Figure 4. Residual spectra (top panels) and spectra in the vicinity of the fc/o, k = 1,...7 pulsation frequency components (bottom 
panels) are drawn. In these panels ~ 2 cd^^ frequency ranges are enlarged. On the X axis ticks are at 2/m separation, corresponds 
to the positions of fc/g. For the sake of lucidity a grid with i/m spacing is also drawn. In the left panels residuals after prewhitening 
with the triplet frequency solution (solution A) are shown. The highest peaks in the spectrum appear at fc/o + 2/„i frequencies. The 
fc/o — 2/m frequencies do not evidently show up in these spectra but in later steps of the prewhitening some of these components also 
emerge with S/N ratio larger than 3. Residuals after prewhitening with the quintuplet frequency solution (solution B) are shown in the 
middle panels. The highest peaks appear at fc/o it 4/m and fc/o ± 12.5/m frequencies. These two series of modulation components are 
unfortunately seriously biased as they are very close to the ±1 cd"'^ alias components of each other. However, after prewhitening the 
data also with 5 frequencies of the fc/o ± 4/m components (septuplet: solution C) as shown in the right panels, it becomes evident that 
the fe/o it 12.5/m components still remain in the spectrum. We conclude therefore, that they are real frequencies, indeed, and not alias 
artifacts of the data sampling. There are other well defined frequency series appearing with common frequency separations, too. One 
is very close to the pulsation components at /m = 0.0019 cd"'^ separation. Four components of this series can be identified: /o — /m, 
2/o ~ /m' ■^fo ~ /m ^'^'i /m- ^^ light curve solution D the septuplet and the fc/o it 12.5/m modulation and also modulation components 
with /m = 0.00197 cd~^ separation are removed. The residuals of light curve solution D are shown in Fig. [5] 



ically the displacements of the other modulation frequency 
components with similar amplitude (see data in Table [5]). 

Some of the other peaks in the residual of the septu- 
plet solution also form series with common frequency sepa- 
ration from the pulsation components. There are remaining 
peaks very close to the pulsation components at /o — /m, 
2/o - /m, 3/o - /m, and /m frequencies with /m = 0.00197 
cd~^. The periods and amplitudes of these frequency compo- 
nents are somewhat uncertain as this modulation period is 
~ 500 days, hardly shorter than the total length of the obser- 
vations. These components can be identified either with an 
additional long period modulation or with residuals caused 
by slight changes of the pulsation period during the obser- 
vations. 

In light curve solution D, besides the quintuplet fre- 
quencies, the fc/o ± 12.5/m and the /m modulation compo- 
nents are also fitted and removed. In Fig. [5] the residual 
spectra of the B, V , and Ic observations are shown after 



prewhitening with light curve solution D. These residual 
spectra are characterized with a broad band low frequency 
signal centred on about the pulsation frequency and with se- 
ries of further modulations at e.g., — 3/m, 5/m, — 11.5/m, and 
13.5/m frequencies. There are also further frequency peaks 
in the ±0.15 cd~^ vicinity of the /o pulsation frequency. 
Though the S/N ratio of the amplitudes of many of these 
frequency components are larger than 3, we stop frequency 
identification at this level as the addition of these frequency 
components does not improve the light curve solution sig- 
nificantly. 

In Table [5] details of light curve solution D are given 
for the B, V, and Ic data. Col. 1 gives the identification of 
the frequencies. Col. 2 lists the frequencies of the 'locked' 
frequency solution (exact linear combinations of the pulsa- 
tion and modulation components). The next 3 columns give 
the differences in the frequency values (A/) of the 'let free' 
solutions for the B, V, and Ic light curves. In the let free 
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Figure 5. The top panels show the residual spectra of the B, V and Iq data prewhitened with frequency solution D (septuplet frequencies 
and modulation components with 12.5/m and f'^ separations). The spectra are zoomed in the vicinity of the fc/o frequencies fc = 1, ...7 in 
the bottom panels. Not regarding the decreasing mean level of the spectra towards longer wavelengths the residuals are similar, peaks at 
— 3/m, 5/m, — 11-5/m, and 13.5/m can be identified. Further significant peaks appear in the |A/| < 3/m vicinity of /q. Though it is very 
probable that many of these frequencies are real signals we stop with frequency identifications here as the addition of these components 
does not significantly improve the light curve solution. 



solutions only the harmonic components of the pulsation 
are at locked frequency values but the best frequency val- 
ues of /o and all the modulation frequency components are 
searched in a nonlinear process. For comparison purposes, 
A//a"/ are given in Cols 6-8 in the three bands. The aj 
error estimates o f the frequencies are calculat e d usi ng the 
formula given by JMontgomerv fc O'Donoghud (|l999l ). It is 
important to note here that, in the case of correlated noise, 
these errors underestimate the true uncertainties of the fre- 
quencies significantly. For most of the identified modulation 
frequency components the frequency displacements are 0.5 
- 4.0 times the calculated (jf values. 

Frequency components with larger displacements 
{Af/af > 3) in each band are: 
/m {Vamp = 0.014); Af{V) = 0.00011 cd-^ 
/o + 2/m (Vamp = 0.004); Af{V) = -0.00055 cd"^ 
4/o - 2/^ {Vamp = 0.003); Af{V) = -0.00064 cd-^ 
/o - 4/m {Vamp = 0.003); Af{V) = 0.00041 cd-^ 
3/o - 4/m {Vamp = 0.002); A/(V') = 0.00090 cd^^ 
fL {Vamp = 0.002); A/(V) = 0.00196 cd^V 

The larger uncertainty of the f^ component arises sim- 
ply from the fact that this modulation period is hardly 
shorter than the time span of the observations. The larger 
displacements of the 4/m components are most probably due 
to their strong ±1 cd~^ alias connections with the frequen- 



cies of the 12.5/m modulation series. The reason why the 
frequency displacement of the other three modulation com- 
ponents are unexpectedly large is unknown. Based on the 
data given in Table [S] we think that there is no serious rea- 
son to assume that the modulation frequencies are not, in 
fact, at their 'locked' positions. 

In Cols 9-14 of Table [5]the amplitudes and phases of the 
B,V , and Ic light curves according to sine term decomposi- 
tion are given. Frequencies with S/N (S/N is defined as the 
ratio of the V amplitude and the mean value of the resid- 
ual spectrum of light curve solution D in the vicinity of the 
given frequency) smaller than 3 in each band are denoted 
by asterisks. The inclusion of these frequencies in the light 
curve solution has minimal effect on the results. However, 
as we are focusing on how accurately the light curve of a 
Blazhko variable can be fitted with the mathematical model 
of equidistant modulation side frequencies, we decided to 
include these low S/N signals also in the frequency solution. 

The rms scatter of the residual light curves of light curve 
solution D in the B, V , and /c bands are 0.026, 0.020 and 
0.015 mag, respectively. Comparing these rms values to the 
residuals of the comparison - check stars' light curves, which 
are in the 0.008 - 0.015 mag range it seems that the residual 
scatter of MW Lyrae is significantly larger than expected. 
The V magnitudes of MW Lyr vary between 12.95 and 14.35 
mag, the brightnesses of the check stars CI, C2, and C3 
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Table 5. Fourier parameters of light curve solution D in the B,V, Ic bands 



Frequency ID 


Frequency 


A/(S) 


A/(V) 


A/(/c) 


^f/Tf 


Af/<Tf 


A//<T^ 


MB) 


«(B) 


A{V) 


#{V) 


Mic) 


Hic) 




[cd-l] 




[10^^ cd-1] 




B 


V 


Ic 


[mag] 


[rad] 


[mag] 


[rad] 


[mag] 


[rad] 


/o 


2.514621 


0.14 


0.39 


0.56 


2.0 


4.2 


3.7 


0.4998 


1.963 


0.3741 


1.919 


0.2302 


1.772 


2/0 


5.029242 


- 


- 


- 


- 


- 


- 


0.2063 


6.211 


0.1585 


6.200 


0.1006 


6.152 


3/0 


7.543863 


- 


- 


- 


- 


- 


- 


0.0968 


4.366 


0.0759 


4.354 


0.0500 


4.335 


4/0 


10.058484 


- 


- 


- 


- 


- 


- 


0.0469 


2.323 


0.0373 


2.334 


0.0251 


2.317 


5/o 


12.573105 


- 


- 


- 


- 


- 


- 


0.0242 


0.305 


0.0199 


0.268 


0.0138 


0.268 


6/o 


15.087726 


- 


- 


- 


- 


- 


- 


0.0165 


4.446 


0.0139 


4.388 


0.0093 


4.389 


7/o 


17.602347 


- 


- 


- 


- 


- 


- 


0.0117 


2.360 


0.0096 


2.394 


0.0069 


2.332 


8/0 


20.116968 


- 


- 


- 


- 


- 


- 


0.0093 


0.324 


0.0074 


0.294 


0.0049 


0.368 


9/o 


22.631589 


- 


- 


- 


- 


- 


- 


0.0065 


4.670 


0.0055 


4.696 


0.0039 


4.574 


lO/o 


25.146210 


- 


- 


- 


- 


- 


- 


0.0053 


2.735 


0.0039 


2.856 


0.0025 


2.592 


ll/o 


27.660831 


- 


- 


- 


- 


- 


- 


0.0031 


0.974 


0.0028 


0.788 


0.0022 


0.940 


12/0 


30.175452 


- 


- 


- 


- 


- 


- 


0.0024 


5.215 


0.0018 


5.329 


0.0010 


5.384 


/m 


0.060437 


8.87 


10.95 


14.93 


4.9 


4.4 


4.3 


0.0193 


3.628 


0.0140 


3.635 


0.0100 


3.558 


/O + /m 


2.575058 


0.06 


0.26 


-0.04 


0.2 


0.7 


-0.1 


0.1217 


3.316 


0.0900 


3.332 


0.0559 


3.352 


2/0 + /m 


5.089679 


0.83 


0.27 


0.24 


2.6 


0.6 


0.3 


0.1067 


1.402 


0.0805 


1.436 


0.0506 


1.501 


3/0 + /m 


7.604300 


0.73 


0.57 


0.83 


1.4 


0.8 


0.8 


0.0648 


6.013 


0.0502 


6.036 


0.0324 


6.103 


4/o + /„ 


10.118921 


3.17 


1.87 


2.58 


4.1 


1.9 


1.7 


0.0451 


4.316 


0.0357 


4.329 


0.0231 


4.352 


5/0 + /m 


12.633542 


2.27 


2.94 


4.63 


1.8 


1.9 


1.9 


0.0274 


2.502 


0.0219 


2.497 


0.0142 


2.491 


6/0 + /m 


15.148163 


9.30 


6.38 


5.81 


4.4 


2.3 


1.4 


0.0165 


0.533 


0.0126 


0.513 


0.0085 


0.525 


7/0 + /m 


17.662784 


6.96 


4.21 


12.06 


2.0 


1.0 


1.9 


0.0098 


4.765 


0.0085 


4.762 


0.0054 


4.763 


8/0 + /m 


20.177405 


11.58 


5.22 


5.85 


2.3 


0.8 


0.6 


0.0069 


2.584 


0.0054 


2.627 


0.0033 


2.676 


9/0 + /m 


22.692026 


6.93 


0.33 


-8.80 


1.2 


0.0 


-0.6 


0.0058 


0.627 


0.0042 


0.595 


0.0024 


0.721 


lO/o + /m 


25.206647 


-9.99 


9.96 


13.16 


-1.2 


1.0 


0.8 


0.0041 


4.919 


0.0035 


4.798 


0.0022 


5.023 


I/O + /m 


27.721268 


-13.31 


19.24 


28.54 


-1.2 


1.7 


1.4 


0.0031 


2.820 


0.0030 


2.894 


0.0016 


2.586 


12/0 + /m 


30.235889 


13.97 


2.98 


65.67 


1.1 


0.2 


1.8 


0.0027 


0.681 


0.0020 


0.987 


0.0010 


0.608 


13/0 + /m 


32.750510 


-9.57 


11.16 


14.54 


-0.6 


0.4 


0.4 


0.0020 


5.425 


0.0014 


5.563 


0.0010 


5.074 


/O - /m 


2.454184 


3.30 


3.74 


2.69 


6.2 


5.2 


2.3 


0.0650 


4.770 


0.0486 


4.787 


0.0299 


4.835 


2/0 - /m 


4.968805 


0.33 


2.23 


0.52 


0.6 


2.9 


0.4 


0.0599 


2.731 


0.0451 


2.775 


0.0278 


2.864 


3/0 - /m 


7.483426 


1.87 


3.54 


3.48 


2.3 


3.4 


2.0 


0.0426 


1.094 


0.0331 


1.109 


0.0204 


1.161 


4/0 - /m 


9.998047 


-1.28 


0.38 


-1.02 


-1.1 


0.2 


-0.4 


0.0287 


5.476 


0.0225 


5.534 


0.0144 


5.575 


5/0 - /m 


12.512668 


4.19 


0.38 


-1.26 


2.2 


0.2 


-0.3 


0.0184 


3.804 


0.0141 


3.763 


0.0094 


3.801 


6/0 - /m 


15.027289 


-4.38 


2.83 


-6.53 


-1.7 


0.8 


-1.2 


0.0132 


1.829 


0.0094 


1.770 


0.0066 


1.768 


7/0 - /m 


17.541910 


6.88 


-1.57 


3.95 


1.6 


-0.3 


0.5 


0.0081 


6.253 


0.0070 


6.239 


0.0047 


6.097 


8/0 - /m 


20.056531 


-19.32 


-14.07 


-3.38 


-2.8 


-1.7 


-0.3 


0.0049 


4.436 


0.0042 


4.422 


0.0026 


4.392 


9/0 - /m 


22.571152 


-31.04 


-21.41 


-5.76 


-2.5 


-1.5 


-0.3 


0.0028 


2.505 


0.0025 


2.473 


0.0017 


2.480 


lO/o - /m 


25.085773 


-20.99 


-1.26 


17.14 


-1.6 


-0.1 


0.7 


0.0026 


0.317 


0.0021 


0.344 


0.0015 


6.220 


2/m 


0.120874 


-2.57 


9.02 


2.30 


-0.7 


2.0 


0.3 


0.0098 


0.335 


0.0075 


0.281 


0.0046 


0.314 


/O + 2/m 


2.635495 


-67.85 


-55.18 


-39.37 


-8.2 


-6.5 


-3.9 


0.0042 


3.216 


0.0041 


3.317 


0.0035 


3.499 


2/o + 2/„ 


5.150116 


-8.37 


-7.21 


7.53 


-1.0 


-0.6 


0.5 


0.0040 


2.666 


0.0030 


2.502 


0.0021 


2.788 


3/o + 2/„ 


7.664737 


9.36 


7.87 


8.36 


2.6 


1.7 


1.4 


0.0098 


1.678 


0.0077 


1.551 


0.0057 


1.512 


4/0 + 2/,„ 


10.179358 


-15.74 


-7.40 


-5.63 


-4.3 


-1.6 


-0.7 


0.0096 


5.573 


0.0075 


5.564 


0.0043 


5.694 


5/0 + 2/„ 


12.693979 


-4.61 


2.84 


0.86 


-1.3 


0.6 


0.1 


0.0095 


4.162 


0.0074 


4.114 


0.0049 


4.172 


6/o + 2/„ 


15.208600 


-23.02 


-17.53 


2.52 


-4.2 


-2.4 


0.2 


0.0063 


2.356 


0.0048 


2.413 


0.0033 


2.384 


7/o + 2/„ 


17.723221 


-7.92 


-0.82 


15.79 


-1.1 


-0.1 


1.2 


0.0048 


0.682 


0.0042 


0.605 


0.0025 


0.709 


8/o + 2/„ 


20.237842 


-16.68 


-18.22 


-16.96 


-1.6 


-1.4 


-0.7 


0.0033 


5.142 


0.0026 


4.945 


0.0013 


4.876 


9/o + 2/„ 


22.752463 


15.03 


17.07 


18.32 


1.0 


0.8 


0.5 


0.0022 


3.303 


0.0016 


3.363 


0.0010 


2.844 


/o - 2/„ 


2.393747 


-2.16 


2.55 


9.01 


-0.3 


0.3 


0.8 


0.0054 


5.786 


0.0044 


5.790 


0.0030 


5.570 


2/o - 2/„ 


4.908368 


-15.89 


-4.45 


-17.27 


-3.5 


-0.7 


-1.6 


0.0077 


5.134 


0.0053 


5.054 


0.0033 


5.351 


*3/o - 2/„ 


7.422989 


38.74 


37.12 


78.85 


2.9 


1.8 


2.3 


0.0026 


3.328 


0.0017 


3.572 


0.0010 


3.798 


4/o - 2/„ 


9.937610 


-67.38 


-63.91 


-68.81 


-7.3 


-5.4 


-3.4 


0.0038 


1.956 


0.0029 


1.977 


0.0017 


2.160 


*5/o - 2/„ 


12.452231 


-49.39 


-56.33 


-70.03 


-3.1 


-2.7 


-2.7 


0.0022 


0.680 


0.0016 


0.445 


0.0013 


0.793 


*6/o - 2/„ 


14.966852 


-112.05 


-123.78 


-101.45 


-4.4 


-4.4 


-2.2 


0.0014 


4.630 


0.0012 


4.592 


0.0008 


4.714 


*7/o -2/„ 


17.481473 


-88.32 


-94.99 


-28.84 


-3.5 


-3.0 


-0.7 


0.0014 


3.053 


0.0011 


3.028 


0.0008 


2.882 


/o + 4/„ 


2.756369 


8.49 


33.32 


8.55 


1.8 


5.3 


1.0 


0.0073 


5.385 


0.0055 


5.408 


0.0040 


5.620 


2/o + 4/„ 


5.270990 


10.58 


32.28 


23.18 


1.2 


2.9 


1.5 


0.0040 


3.604 


0.0031 


3.664 


0.0022 


3.654 


/o - 4/„ 


2.272873 


98.12 


41.49 


70.76 


7.6 


3.7 


3.6 


0.0027 


1.948 


0.0031 


2.402 


0.0018 


2.443 


2/o - 4/„ 


4.787494 


50.96 


54.76 


-3.99 


6.0 


3.8 


-0.2 


0.0041 


5.921 


0.0024 


5.884 


0.0021 


6.080 


*3/o -4/„ 


7.302115 


88.66 


90.46 


97.93 


7.5 


5.1 


4.0 


0.0029 


5.100 


0.0020 


4.930 


0.0014 


5.341 


/o + 12.5/,„ 


3.270084 


-45.30 


-20.18 


-12.27 


-3.7 


-1.7 


-0.8 


0.0028 


3.050 


0.0030 


2.960 


0.0024 


2.956 


2/o - 12.5/,„ 


4.273780 


9.48 


39.58 


50.07 


1.2 


3.7 


2.4 


0.0044 


3.821 


0.0033 


3.783 


0.0017 


4.018 


*2/o + 12.5/,„ 


5.784705 


-143.83 


-36.28 


-101.75 


-9.4 


-2.8 


-5.3 


0.0023 


1.096 


0.0026 


1.336 


0.0018 


1.599 


3/o - 12.5/,„ 


6.788401 


-17.73 


-74.47 


-28.12 


-2.6 


-9.0 


-2.1 


0.0051 


2.557 


0.0042 


2.473 


0.0025 


2.346 


fo-fL 


2.512651 


30.64 


39.23 


25.88 


7.1 


5.2 


2.4 


0.0081 


2.906 


0.0046 


2.718 


0.0033 


2.567 


2/o - /,;. 


5.027272 


29.86 


-12.99 


-99.47 


4.4 


-1.5 


-5.9 


0.0051 


0.709 


0.0039 


0.484 


0.0020 


0.600 


3/o - fL 


7.541893 


7.10 


-4.20 


-48.85 


0.8 


-0.5 


-3.2 


0.0038 


5.460 


0.0039 


5.365 


0.0023 


5.214 


fL 


0.001970 


110.33 


196.07 


195.70 


11.4 


12.4 


9.3 


0.0036 


0.717 


0.0022 


0.539 


0.0016 


0.127 



nearly equal to the variables maximum, mean and minimum 
brightnesses, respectively (see the magnitudes and colours 
of the comparison and check stars in Table [TJ . The mean 
A{B - V) between MW Lyr and CO is -0.25 mag, while 
the A{B — V) colour differences of the three check stars are 
between 0.02 and 0.28 mag. Neither the brightness nor the 



colour differences can account for the larger scatter of the 
residual light curve of MW Lyr. 

Fig. |6] shows some examples how hght curve solutions 
A, B, and D fit the observations on different nights. It can 
be seen that the observations deviate systematically even 
from the fit of solution D that involves 66 frequency com- 
ponents. The deviations are systematic, they are centred on 
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V residual spectrum of MW Lyr - CO: light curve solution D removed 



Figure 6. V light curves and the triplet (dotted blue line), quin- 
tuplet (dashed black line), and light curve solution D (solid red 
line) fits to the data are shown on some representative nights. On 
the X axes tick marks are at 0.005 d on each plot. Bellow the 
light curves the residuals if light curve solution D is removed are 
shown. The residuals show smooth light curves with systematic 
deviations from zero. Large residuals appear at around minimum 
and maximum light and on the rising branch. However, as e.g., on 
JD 2454200 the residual is continuously small but negative on the 
upper part of the descending branch. These systematic deviations 
explain the large rms. scatter of the residual light curve. 



the minimum-rising branch-maximum phase of the pulsation 
with different shapes. The residual mimics the behaviour of 
the large amplitude modulation on a much smaller scale, but 
without any definite periodicity. 

Just for a trial we have included further 30 frequencies 
in the fit (consecutive prewhitening with the highest peaks 
appearing in the spectra shown in Fig.O hght curve solution 
E) to see whether these systematic deviations can or cannot 
be eliminated with further distinct frequency components. 
Though the rms of e.g., the V hght curve has been reduced 
to 0.018 mag this way, the systematic deviations of the hght 
curves shown in Fig. |6] have hardly decreased. 

Fig.[7]compares the residual spectra of the V hght curve 
prewhitened with hght curve solution D and E with the spec- 
tra of the V hght curves of the CI, C2 and C3 check stars. 
The differences are striking. The mean level of the residual 
spectrum of hght curve solution E is still 2 — 3 times larger 
than that of the spectra of the check stars. 

Most probably stochastic and/or chaotic behaviour of 
the modulation bring forth the enhanced residual signals. It 
seems that the hght curves of large modulation amplitude 
Blazhko variables cannot be modelled with the required ac- 
curacy with the Fourier sum of finite number of frequency 
components. 



3.2 The maximum brightness — maximum phase 
variation 

The modulation of the pulsation hght curve of MW Lyr can 
be also followed using the maximum brightness, maximum 
phase data given in Table [3] The maximum brightness and 
maximum phase values of the V hght curve are shown in 
the left panels of Fig. [S] folded with the modulation period. 
None of these plots can be fitted with a single sine wave. 
The maximum phase data show a highly asymmetric shape, 
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Figure 7. Comparison of the residual spectra of the V light curve 
of MW Lyr with the spectra of the check stars' light curves. 
The identifications and magnitudes of the CO comparison and 
CI, C2, C3 check stars arc given in Table [T] Two residual spec- 
tra of MW Lyrae arc shown. After removing the 66 frequencies 
of light curve solution D significant peaks in the spectrum still 
appear (see also Fig. [Sj- Removing further 30 frequencies from 
the data the residual spectrum is still at 2-3 times higher level 
than the noise spectrum of the check stars. As the brightness of 
MW Lyr varies between 13.10 mag and 14.30 mag, while CI, C2, 
and C3 have brightnesses of 12.99, 13.92 and 14.36 mag, respec- 
tively, brightness differences cannot account for the large residual 
signal of MW Lyrae. 



while there is a bump on the rising branch of the maximum 
brightness data. For an accurate fit, both plots need at least 
6th order Fourier sums. Table |6] lists the amplitudes of the 
6th order Fourier components of the maximum brightness 
and maximum phase fits. The amplitudes change with in- 
creasing order, they follow a similar trend for the maximum 
brightness data as the amplitudes of the sidelobe frequen- 
cies of the Fourier spectrum of the hght curve do: the 1st, 
2nd and 4th order components have pronounced amplitudes 
while the amplitudes of the 3rd, 5th and 6th order compo- 
nents are small. On the contrary, the amplitudes of the 3rd, 
4th and 5th order components of the maximum phase fit are 
very similar. The 0.022 mag and 0.0066 phase (0.0026 days) 
rms scatter of the residuals of the maximum brightness and 
maximum phase fits are, however, too large as we estimate 
that the data are accurate within ±0.02 mag and ±0.0025 
d, respectively. 

The residual spectra of the maximum brightness and 
maximum phase data prewhitened with the 6th order fits 
of the modulation frequency are shown in the right panels 
of Fig. [8l The residual spectrum of the maximum brightness 
data shows a peak with 0.015 mag amplitude at 0.7559 cd~^. 
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Figure 9. V maximuin brightness vs. maximuin phase data of the 88 observed maxima. Data for the consecutive 50-70 days segments 
(3-5 Blazhko cycles) are shown by different symbols in the left panel. No systematic change with time in the Vmax — Phascmax plot 
is evident. For comparison, in the I, and II panels synthetic V^max — Phasemux curves according to the light curve solution A (triplet) 
and D (that involves the septuplet solution of the modulation plus 4 modulation components with 12.5 /m separation, and 4 modulation 
components with f^ = 0.0019 cd~^ modulation frequency) are plotted. In panel III the fitted curve corresponds to the harmonic fits 
to the maximum brightness and maximum phase data with the frequencies given in Table |6] The observed V^ax — Phassmax data are 
overplotted on the synthetic solutions, the open symbols denote less accurate data. The scatter of the observations is much larger than 
observational uncertainties can explain as data are accurate within ±0.02 mag and ±0.003 d ±0.0075 phase ranges. 
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Figure 8. Maximum V brightness and maximum phase values vs. 
Blazhko phase are plotted in the top-left and bottom-left panels, 
respectively. Sixth order Fourier fits to the data are also drawn 
in these plots. The data are prewhitened with these fits. The 
right-top panel shows the amplitude spectrum of the residual 
of the maximum brightness data. The highest peak appears at 
/ = 0.7559 cd~^, which equals to 12.5/m within the uncertainty. 
Prewhitening the data also with this frequency, there is no further 
significant peak in the residual spectrum (middle-right panel). 
The residual spectrum of the prewhitened maximum phase data 
(bottom-right panel) does not show the / = 0.7557 cd^^ sig- 
nal. There is no peak appearing at the same frequency in the 
maximum phase residual (prewhitened with the modulation fre- 
quency and its harmonics) and the maximum brightness residual 
(prewhitened with the modulation frequency and its harmonics 
plus 12.5 /m) spectra. 



This frequency equals within the uncertainty to the 0.7555 
cd~^ frequency value of 12.5/m. Though this modulation 
component itself does not appear in the spectrum of the 
Ught curve, it has measurable amplitudes at the fc/o± 12.5/m 



sidelobe positions as discussed in the previous Section. There 
is no sign of this modulation component in the maximum 
phase data. Accordingly, the modulation connected to the 
12.5/m frequency is dominantly amplitude modulation. 

The removal of the 12.5/m component from the max- 
imum brightness data lowers the residual scatter to 0.020 
mag only, which is still too high to be explained with ob- 
servational inaccuracy. Both the maximum brightness and 
maximum phase data reflect the imperfection of modelling 
the modulation using the Fourier sum of discrete frequency 
components. Most probably the observations cannot be 
traced with the expected accuracy using finite number of 
strictly periodic signals. 

Plotting the maximum brightness vs. maximum phase 
data, as Fig. [5] shows, significant scatter appears around 
a triangular shape curve. No evolution of the data with 
time account for the scatter, data from different segments of 
the observations are equally scattered without any system- 
atics. Synthetic maximum brightness and maximum phase 
curves are drawn according to the triplet light curve solu- 
tion, light curve solution D, and maximum brightness and 
maximum phase fits involving 6 harmonics of the modula- 
tion and 12.5/m for the maximum brightness data in the 
panels I, II and III in Fig. [^ respectively. Many observed 
maximum data are out of the ranges of any of the model 
fits, that cannot be explained by data inaccuracy. 



3.3 The mean light curve 

Interesting questions arise in connection with the mean light 
curves of variables showing light curve modulation. Namely, 
we still do not know whether 

a) the mean light curve is or is not the same as the light 
curve of a star with the same physical parameters but not 
showing Blazhko modulation, 

b) the actual shape of the light curve in any phase of the 
modulation corresponds to the mean light curve. 

Our extended and accurate data make it possible to 
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Table 6. Fourier amplitudes of the fits to the maximum bright- 
ness and maximum phase data 



frequency 



maximum brightness 
Amp [mag] 



maximum phase 
Amp [pulsation phase] 
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0.0047 
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Figure 10. Comparison of the mean light curves defined as a) 
the Fourier fit of the pulsation frequency and its harmonics to the 
y observations from the first season (solid red line) , b) the Fourier 
parameters of the pulsation components are taken from a triplet 
solution of the data of the first season (dashed black line) , c) the 
Fourier parameters of the pulsation components are taken from 
light curve solution D of the entire dataset (blue dot-dashed line). 
The residual spectra in the vicinity of the pulsation frequency 
and its first harmonic component are shown in the right panels. 
In method a the modulation side components are also removed 
but in a second step after the pulsation frequencies have been 
already removed from the data. Note the smaller amplitude of 
the fitted light curve and the high residual signals at /q and 2/o 
when applying method a. 



find the correct answers to these questions. However, first 
we have to define the mean pulsation light curve of Blazhko 
variables correctly. 

The simplest way to define the mean light curve is to 
fit all the data with a high enough order Fourier sum of the 
pulsation frequency and its harmonics (method a). How- 
ever, this procedure may yield a mean light curve that is 
biased by the uneven data sampling, even in the case of an 
extended data set. The typical procedure of light curve anal- 
ysis of Blazhko stars is prewhitening first with the pulsation 
frequency and its harmonics (with the mean light curve de- 
fined above), and then, the prewhitened data are analysed 
to identify the modulation components. As a result of the 
incorrect shape of the removed mean light curve during this 
procedure false signals in the vicinity of the pulsation com- 
ponents emerge. A better approximation of the mean light 
curve is gained by a fit that takes the pulsation and mod- 
ulation frequency components simultaneously into account. 
The mean pulsation light curve according to a triplet solu- 



tion light curve fit (method 6) may significantly differ from 
the mean light curve defined by method a. 

Fig. [TD] demonstrates the differences between the 
fits and the residual spectra if the data are fitted and 
prewhitened by the pulsation and modulation components 
in consequtive steps and simultaneously. In this Figure data 
from the first season of the observations are shown. This 
data set, which contains 3900 data points from 120 nights is, 
however, more dense and extended than any previous pho- 
tometric observation of a Blazhko variable. Though these 
data cover the whole pulsation light curve in each 0.05 phase 
of the modulation, the differences between the mean light 
curves defined by method a, and by method 6 are signifi- 
cant. 

The fitted mean light curves shown in Fig 1101 are: 

(i) the Fourier fit to the data with the pulsation frequency 
and its harmonics {method a); 

(ii) the fit taken from the simultaneous triplet frequency 
solution of the data {method b); 

(iii) the pulsation light curve taken from the complete 
light curve solution of the entire data set (light curve solu- 
tion D as given in Table [SJ. 

Method a gives substantially smaller amplitude around 
minimum and maximum light due to the oversampling of 
the small amplitude phase of the modulation in the data, 
than the fits either from method b or from the light curve 
solution D. The latter fits are very similar, showing that with 
the inclusion of the most prominent modulation components 
(triplets) in the light curve solution quite a reliable mean 
light curve can be gained even if the data sampling is biased. 

The right panels in Fig. [10] show the residual spectra of 
the three fits in the vicinity of the pulsation frequency and 
its first harmonic. In order to make the residuals compara- 
ble, the modulation side frequencies are also removed in a 
second step in method a. These residual spectra show high 
amplitude signals at /o and 2/o, which also points to the 
inadequacy of prewhitening the data in consequtive steps. 

We thus conclude that one has to be cautious how to 
define the mean light curve of a Blazhko variable as from 
method a and method b different results emerge. In real ob- 
servations data sampling is always somewhat unevenly dis- 
tributed. As a consequence, the mean light curve defined 
without taking the modulation components also into ac- 
count, may give an incorrect result for variables showing 
any type of light curve modulation. 

From here on, we use the Fourier parameters of the 
pulsation components given in Table \S\ to define the mean 
pulsation light curve. 

3.4 Light curve changes during the Blazhko cycle 

Utilizing the full coverage of the pulsation period in each 
phase of the modulation in our data set we can reliably 
compare the mean light curve to the light curves in dif- 
ferent phases of the Blazhko modulation. In Fig. [11] the V 
light curves of MW Lyr are shown for 20 bins of the mod- 
ulation cycle. The scatter of these light curves can be par- 
tially explained by the regular light curve changes that take 
place even in 0.05 phase intervals of the modulation. Modu- 
lation frequencies that are not integer multiplets of /m (e.g., 
12. 5/m and f^) also result in enhanced scatter of the data 
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Table 7. Fourier parameters of the V light curves in the 0.05 phase bins of the modulation 



Bl phase Order 



<V> 



[rad] 



$41 



A(/o) A(2/o) 



A(3/o) 
[mag] 



A(4/o) A{5/o) 



0.00-^ 


-0.05 


14 


0.015 


2.031 


228 


1.691 


2.423 


5.173 


1.483 


4.146 


0.459 


0.244 


0.148 


0.092 


0.059 


0.05-^ 


-0.10 


14 


0.024 


2.024 


277 


1.737 


2.430 


5.136 


1.465 


4.166 


0.460 


0.247 


0.142 


0.087 


0.052 


0.10- 


-0.15 


13 


0.032 


2.017 
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1.792 


2.417 


5.095 


1.364 
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0.457 


0.234 
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0.077 
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0.15- 


-0.20 


13 


0.020 
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0.209 
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13 
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5.030 
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-0.30 


13 
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2.053 
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Figure 11. V light curves of MW Lyrac in 20 different phases of the Blazhko cycle are shown. The folded light curve of the complete 
data set is shown in gray colour. 
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Figure 13. Left panel: Fitted curves of tlie V light curve of 
MW Lyr in 20 bins of the Blazhko cycle. Each light curve is phase 
shifted with the 'J'(/o) phase value. Right panel: Time transformed 
V light curve of MW Lyrae. The time transformation was defined 
from a second order harmonic fit to the "l?(/o) phases in different 
phases of the modulation. 



Figure 12. Fourier parameters of the V light curves in 20 phase 
bins of the Blazhko cycle are plotted. The open circle in the $51 
plot denotes a very uncertain data point, the amplitude of the 
5/0 component is only 0.002 mag in this phase bin. The ampli- 
tudes and the phases are plotted on the same scales. Note that 
the Fourier phases have opposite sign as the directly measured 
phase of the maximum brightness shown in Fig |8] The ampli- 
tudes of the fc/o {k = 1,...5) components and the 3>(/o) phase 
show smooth, sinusoidal variations with the smallest phase value 
when the amplitudes are the highest and with the largest phase 
value when the amplitudes are the smallest. On the contrary the 
epoch independent phase differences ('I'fci) show surprisingly dif- 
ferent behaviour. There are small if any changes in 'I'2i> and the 
variation in <3?3i is also very small taking into account the signifi- 
cant changes in the light curves' shapes as shown in Fig. 1111 <I'4i 
and <&5i vary during the Blazhko cycle showing complex changes 
on the time scale of about half the modulation period. For com- 
parison, Fourier parameters of the mean light curve as given in 
Table \5\ are also drawn by dashed lines in the plots. While the 
<t>(/o) value of the mean light is close to the average of their 
observed values in the different Blazhko phases, the phase differ- 
ences in each phase of the modulation are larger than the phase 
differences of the mean light curve. The Fourier amplitudes of the 
mean light curve are 0.01 — 0.02 mag fainter than the averages of 
their values in different phases of the Blazhko cycle. 



phased with the modulation period. Moreover, our experi- 
ence, that the hght curve cannot be fitted with the required 
accuracy supposing regular modulations with different mod- 
ulation periods as discussed in Sect 3.1 and 3.2, means that 
this irregular character of the modulation adds some extra 
noise to the hght curves in the different phase bins, as well. 

Notwithstanding these effects, the hght curves in the 
different phase bins are well defined, and can be character- 
ized by the Fourier amplitudes and phases of the pulsation 
frequency and its harmonics. 

The changes in the Fourier amplitudes and phases 
(phase differences) of /o and its lower harmonics and also the 
average values of these parameters are listed in Table [7] For 
comparison, the last line gives the corresponding parameters 
of the mean V hght curve. The first five columns in Table [7] 
give the phase bin, the order of the Fourier sum fitted to the 
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Figure 14. Phase relation of the amplitude (dots) and period 
(red line) changes of MW Lyrae during the Blazhko cycle. The 
amplitude and period variations are measured as the amplitude 
and as the derivative of the phase variation of the /o pulsation 
frequency component, respectively. 



data, the residual scatter of the fit, the intensity weighted 
mean magnitude, and the number of data points belonging 
to the given bin, respectively. The phases of the /o pulsa- 
tion frequency can be read from the sixth column (initial 
epoch is 2 453 887.00) while the next four columns list the 
epoch independent phase differences $^1 (fc = 2, ...5). The 
amplitudes of the fc/o {k — 1, ...5) components are given in 
the last five columns. The residual scatters of the fitted har- 
monic functions to the data in the different phase bins of 
the modulation are within the 0.015-0.032 mag range. The 
average of the rms values of the fits is 0.022 mag, slightly 
larger than the rms of the hght curve solution D of the entire 
data set. 

The data fisted in Table [7] and plotted in Fig. [12] show 
that there are major differences between the Fourier param- 
eters of the mean light curve (the Fourier parameters of the 
fc/o components of the full hght curve solution given in Ta- 
ble[S} and the average values of the same Fourier parameters 
of the fight curve fits in different phases of the Blazhko cy- 
cle. The amplitudes of the fc/o components of the mean hght 
curve are systematically 0.01 — 0.02 mag smaller than the 
averages derived from the fits of the individual light curves. 
The $21, "i&si, $41 and $51 phase differences of the mean 
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light curve solution are 0.16, 0.17, 0.35 and 0.78 rad smaller 
than the averages of the corresponding parameters of the 
light curves in different phases of the Blazhko cycle. 

These large differences make it unambiguous that in 
each phase of the Blazhko modulation the light curve of 
MW Lyr differs from its mean light curve. 

The $fei phase differences show quite a surprising be- 
haviour. In spite of the large amplitude of the Blazhko mod- 
ulation of MW Lyrae, $21, and $31 hardly vary, while the 
changes in $41 and <l>5i show complex behaviour with cy- 
cle length half of the modulation cycle. For the higher order 
components the amplitudes at around Blazhko minimum are 
so small that the errors in the phases become too large to 
make any firm conclusion about their variations. 

The light curve changes during the Blazhko cycle can 
be summarised as follows: 

Light curve changes connected to the amplitude variations: 
The amplitudes of the fc/o frequencies show parallel 
changes but the amplitudes of the higher order compo- 
nents decrease more drastically than the amplitudes of the 
lower order components. For example, the amplitude ratio 
A{fo)/A[f^) is about 7 in the three largest amplitude phase 
bins, but this ratio is as large as 50-150 for the three smallest 
amplitude phase bins. 

The light curves at around maximum amplitude phase 
of the modulation can be fitted accurately with 10-15 har- 
monic components of /o. The light curves are much more 
sinusoidal in the small amplitude phases, they can be fitted 
with 6-9 harmonic components with the required accuracy. 

Light curve changes connected to the phase variations: 

The phases of the lower order harmonic components 
show very harmonized changes as indicated by the small 
variations of the lower order phase differences. This means 
that the changes in the light curves' phase can be charac- 
terized basically by one parameter, with the phase of the /o 
pulsation frequency, $(/o). 

The $41 and $51 phase differences show double wave 
curves indicating that the time scale of their variation is 
about half of the period of the modulation. The appearance 
of the 2/ni, 4/m frequency components in the Fourier spec- 
trum of the whole data set is probably connected with this 
double periodic behaviour of the higher order phase differ- 
ences. 

We have also tested how accurately the observations 
can be fitted with a simple mathematical model which de- 
scribes the amplitude and phase changes of the different 
order pulsation components with different order Fourier se- 
ries. We have found that such a light curve solution fails 
to fit the observaions with similar accuracy as the Fourier 
sum of the pulsation and modulation side lobe frequencies 
involving even a smaller number of parameters. This result 
means that the modulation has a very complex behaviour 
that could not be described by a mathematical model of 
amplitude and phase modulations of harmonic functions. 

Phenomenologically, amplitude modulation occurs if 
there are observed changes in the brightnesses of the light 
curve maxima, while phase modulation manifests itself as 
a missing fix point on the rising branch of the folded light 
curve. 

Exploiting the slight changes in the lower order phase 
differences, we can 'harmonize' the phases of the light curves 



with a simple phase correction. Shifting the fitted light 
curves of the different bins by their $(/o) values a surpris- 
ingly coherent light curve series emerge as shown in Fig. 1131 
The pronounced fix point occurring on the rising branch 
in this figure validates our simple treatment separating the 
phase modulation of the light curve from the amplitude 
modulation by correcting the phases with the $(/o) values. 

A similar procedure can be applied on the whole data 
set, as well. Fitting the <i>(/o) values of the 20 phase bins with 
a 2nd order Fourier sum we can define a continuous function 
of the phase variation. It is supposed that by transforming 
the times of the observations according to this function we 
'get rid of the phase modulation component of the modu- 
lation. As the right panel in Fig. [13] shows, this is indeed 
the case, the light curve of the time transformed data shows 
very regular amplitude modulation with a pronounced fix 
point on its rising branch. The phase difference of the am- 
plitude peaks of the time transformed data is much smaller 
than in the original data. Its range is consistent with the 
value expected from the amplitude variation. 

Without starting the light curve analysis from the be- 
ginning using this phase corrected data set, we have only 
checked 

• whether or not the same frequency components occur in 
the spectrum of the time transformed data as in the original 
data set, 

and 

• how the rms scatter of the time transformed data com- 
pares to the rms residual of the original data. 

All the frequencies but 2/o — 4/m and 3/o — 4/m listed 
in Table [S] can be identified in the residual spectrum of the 
time transformed data. The fm and 2/m modulation side- 
lobe components can be detected up to higher order harmon- 
ics as contrasted with the spectrum of the original data. 

The residual scatter of the original V light curve was 
0.020 mag, the rms scatter of the time transformed data 
remains the same, 0.020 mag if the modulation components 
up to the appropriate order are taken into account. The 
Fourier parameters of the mean pulsation light curve of the 
time transformed data (given in the last line in Table [7| 
equal within their error ranges with the average values of the 
Fourier parameters of the observed light curves in different 
phases of the Blazhko cycle. The residual light curve of the 
time transformed data also shows large deviations at around 
minimum-rising branch-maximum phases of the pulsation. 

Without finding the correct explanation of the Blazhko 
phenomenon we cannot decide which characterization of the 
light curve is correct: the Fourier analysis of the light curve 
as it is, or the separation of the phase and amplitude mod- 
ulation components of the light curve modulation by an ap- 
propriate transformation of the phases (times) of the obser- 
vations. The simplicity and the low number of independent 
parameters involved in the time transformation applied sug- 
gest that this new treatment of the light curve may lead to 
a headway in the study of Blazhko variables. 

We also remark that the time transformation can be ex- 
plained as a continuous change in the pulsation period dur- 
ing the Blazhko cycle as lStotherd (120061 ) has also interpreted. 
The full range of period change determined from the deriva- 
tive of the phase shift curve is 0.006 days, i.e. SP/P = 0.015. 
The period of the pulsation is about 0.401 days around 
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Blazhko phase 0.65 and 0.395 days around Blazhko phase 
0.30 as Fig. [14] shows. The a mphtude o f the period vari- 
ation is somewhat larger than IStothera (|2006h derived for 
RR Lyrae itself, but keeping in mind that for RR Lyr tem- 
poral periods were determined for some days long intervals 
of the observations, most probably a reduced value of its 
real period change was found. 



4 CONCLUSIONS 

The photometric observations of MW Lyr analyzed in this 
paper comprise the most extended and accurate data set 
of a Blazhko variable ever obtained. Utilizing this unique 
opportunity a detailed and circumspect phenomenological 
description of the modulation is given, which may provide 
crucial information to find the correct explanation of the 
phenomenon. 

The main results and conclusions of the analysis are the 
foUowings: 

• In the Fourier spectrum of the light curve besides the 
fc/o ±/m triplet frequencies fc/o ±2/m quintuplet, and fc/o ± 
4/m septuplet components also appear. 

• Both /m and 2/ni frequencies can be detected in the 
spectrum. 

• Frequency components at fc/o ± 12.5/m are detected. If 
these frequencies are not 'just by chance' at ±12.5/m sepa- 
rations but somehow they are indeed connected to the main 
modulation frequency, then it is a great challenge to find an 
answer to their origin. 

• Modulation with 12.5/m frequency can also be detected 
in the maximum brightness data but not in the maxi- 
mum phase observations. Consequently, the modulation con- 
nected to the 12.5/m frequency is dominantly amplitude 
mo dulation. 

• IStotherj (|2006l ) mentioned the lack of large amplitude 
modulation around minimum phase of the pulsation as a fail- 
ure of his model. The amplitude of the modulation around 
minimum phase of the pulsation in MW Lyr is, however, 
commensurable with the amplitude of the modulation in 
maximum brightness. On the contrary, the modulation of 
MW Lyrae is more strictly periodic and regular than it 
would be expected to be if the triggering mechanism behind 
the modulation were the cyclic weakening and strengthen- 
ing of the turb ulent convection in the ionization zones as 
IStothersI (|2006l ') proposes. 

• Though the modulation shows high degree of regularity 
both in the phase (period) and in the amplitude changes, the 
light curve cannot be fitted with the required accuracy even 
with 66 identified and further 30 frequencies appearing in 
the residual spectrum. Significant deviations in the residual 
light curve are concentrated at the minimum-rising branch- 
maximum phase of the pulsation. These residuals, however, 
do not show any periodicity, most probably they can be 
explained with some stochastic and/or chaotic behaviour of 
the modulation itself. 

• The mean pulsation light curve defined by the Fourier 
parameters of the fc/o pulsation frequency components of the 
full light curve solution differs significantly from the light 
curve in any phase of the modulation. Especially the $*,! 
phase differences of the mean light curve are discrepant, they 



are out of the range of the phase difference values measured 
in any phase of the Blazhko cycle. 

• The light curves in the small and large amplitude phases 
of the modulation can be fitted with 6-10 and 11-15 order 
harmonic fits, respectively. If nonlinear effects (e.g., shock 
waves) account for the occurrence of the higher order har- 
monic components of the pulsation, then their diminishing 
amplitudes in the low amplitude phase of the modulation 
may indicate that these nonlinear effects are not so impor- 
tant in_thisphaseof the modulation. 

• IStotherd (|2006r ) proposed that the enhanced convection 
lowers the pulsation amplitude in the small amplitude phase 
of the modulation, while the phase relation betwen the pe- 
riod changes and the amplitude variations depends on the 
physical parameters of the individual variable in a complex 
way. The pulsation period of MW Lyr is the longest about 
1-2 days (0.1 phases of the modulation cycle) later than the 
minimal amplitude phase, while it is the shortest about 3- 
4 days (0.2 phases of the modulation cycle) later than the 
maximal amplitude phase. 

• The small variations in the phase differences of the 
lower order harmonic components of the light curves dur- 
ing the Blazhko cycle indicate that the phase variations of 
these frequency components are very coherent. If the mod- 
ulation were caused by the intera ction of close radial mode 
and nonradial mode frequencies (|Dziembowski fc Mizerskil 
120041 ) ■ then no such phase coherency would be expected. 

• The modulation of the light curve of MW Lyr can 
be separated into amplitude and phase modulation compo- 
nents using only one parameter, the phase of the /q pulsa- 
tion frequency in each phase bins of the modulation. This 
is in agreement with the explanation of the Blazhko effect 
with per i od an d amplitude changes as recently proposed by 
IStothersI (|2006l '). 

It would also be important to know how common the de- 
tected properties of the light curve modulation of MW Lyrae 
are. Only further similar observations of other Blazhko stars 
can give an answer to this question. 
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